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We extend the work of Oppenlreimer & Synder to model the gravitational collapse of a star to a 
black hole by including quantum mechanical effects. We first derive closed-form solutions for classical 
paths followed by a particle on the surface of the collapsing star in Schwarzschild and Kruskal 
coordinates for space-like, time-like and light-like geodesics. We next present an application of these 
paths to model the collapse of ultra-light dark matter particles, which necessitates incorporating 
quantum effects. To do so we treat a particle on the surface of the star as a wavepacket and 
integrate over all possible paths taken by the particle. The waveform is computed in Schwarzschild 
coordinates and found to exhibit an ingoing and an outgoing component, where the former contains 
the probability of collapse, while the latter contains the probability that the star will disperse. These 
calculations pave the way for investigating the possibility of quantum collapse that does not lead to 
black hole formation as well as for exploring the nature of the wavefunction inside r = 2 M. 

PACS numbers: 


A. Introduction 

Black holes play a pivotal role in the evolution of the 
universe providing an important test laboratory for gen¬ 
eral relativity. Within classical general relativity, Oppen- 
heimer & Synder modelled the gravitational collapse of 
star to a black hole by approximating the star with a uni¬ 
form sphere of dust (hereafter O-S model) )lj. This model 
provides an analytic solution for stellar collapse that con¬ 
nects the Schwarzschild exterior of a star to a contracting 
Friedmann-Robertson-Walker (FRW) interior. Once the 
surface has passed within r = 2 M, no internal pressures 
can halt the collapse and all configurations collapse to 
a point-like singularity at r = 0. The general features 
of this toy collapse model have been examined by many 
authors HU. 

In the standard formalism from [3j, the stellar surface 
is considered to be initially at rest. Here we consider 
configurations with all possible initial velocities. We de¬ 
rive closed-form solutions for the equations of motion 
in Schwarzschild and Kruskal coordinates for space-like, 
time-like and light-like geodesics. 

As an example application of our closed-form solutions 
for the classical O-S model with non-zero initial veloci¬ 
ties, we consider the macroscopic collapse of a spheri¬ 
cally symmetric sphere of dust composed of ultra-light 
particles. To approximate quantum effects, the radius 
of the star is approximated by a Gaussian wavepacket 
that is initially centred far from r = 2 M. Its evolu¬ 
tion is then followed via a simple path integral approach 
that extends the results of Redmount and Suen 7] from 
a relativistic free particle to a particle constrained by 
non-trivial gravity. Gravitational collapse incorporating 
approximate treatments of quantum mechanics has been 
considered by a variety of authors EHm. They involve 


non-equivalent ways of quantization that often produce 
physically different results mi. Our path-integral ap¬ 
proach approach is simpler, and can be easily compared 
to the assumption that the particle obeys the relativistic 
Schrodinger equation. 

In the classical model a star is idealised as a collaps¬ 
ing self-gravitating dust sphere of uniform density and 
zero pressure where the constituent particles have the at¬ 
tributes of classical dust: each particle is assumed to be 
infinitesimal in size and to interact only gravitationally 
with other matter. The inclusion of quantum mechani¬ 
cal effects lifts some of these assumptions allowing for the 
possibility that some configurations will not collapse to 
black holes but will disperse or even form stable new con¬ 
figurations. Quantitatively, quantum treatment is neces¬ 
sary in macroscopic stellar collapse when the action S is 
of the order h HE). In our case, the S/h ~ 1 condition 
corresponds to mM ~ Ti. For macroscopic black holes of 
masses M = 1 Mq — 10 9 Mq, this implies that quantum 
treatment is necessary for ultra-light constituent parti¬ 
cles of m ~ 10~ 10 — 10 -19 eV. In nature, such ultralight 
particle could be dark matter. Dark matter comes close 
to the attributes of classical dust, and some dark matter 
clouds may be dense enough to collapse to black holes. 
Particles as light as 10 -22 —10 -23 eV have been proposed 
as constituents of dark matter halos piffl - SHl . 

To approximate quantum effects of the collapse of such 
a halo comprised of ultra-light particles, we start with an 
initial wavefunction that represents the position of the 
particle on the stellar surface, and is at first far from 
2 M. To study the propagation of the wavefunction, we 
integrate over all possible paths taken by the particle. 
At a given time, the outgoing wavefunction comprises 
the probability that the star disperses, and the ingoing 
wavefunction the probability that it collapses. We com¬ 
pute the propagator in analytical form for a particle on 


2 


the surface of the star in the WKB approximation, and 
compare this to the limited assumption that the par¬ 
ticle obeys the relativistic Schrodinger equation. Our 
equations reduce to a free particle case when the mass 
of the star is zero; then the solution to the relativistic 
Schrodinger equation is the exact representation of the 
wave function [7]. In the more general case of a par¬ 
ticle on the surface of a star, this representation is no 
longer correct. However, the comparison is instructive. 
As expected, we observed that the WKB and relativistic 
Schrodinger approximations are out of step at early and 
late times, and appear to converge towards one another 
at intermediate times. 

In Schwarzschilcl coordinates, we can follow the evolu¬ 
tion of the surface of the star only until the formation 
of an apparent horizon due to the r = 2 M coordinate 
singularity. In Kruskal coordinates, one can continue to 
study the evolution of the star inside the apparent hori¬ 
zon (r < 2 M). The paths we have derived here include 
time-reversing space-like paths, which allow for the pos¬ 
sibility of extraction of information from inside the hori¬ 
zon to the outside. These paths turn inside the horizon 
r > 0, and head toward r = 2 M. In future work, explicit 
geodesic equations in Kruskal coordinates may be used 
as a stepping stone to model behaviour inside r = 2 M. 

The rest of the paper is structured as follows: in Sec¬ 
tion I we describe the Oppenheimer-Snyder model. Sec¬ 
tion II includes an overview of the classical action and 
paths, a computation of the Oppenheimer-Snyder limit 
and the derivation of the classical paths in Schwarzschild 
coordinates. Section III details the classical paths in 
Kruskal coordinates, which comprise a starting point for 
future work that explores the collapse inside r = 2 M. As 
an example application of the closed form geodesics, Sec¬ 
tion IV describes the quantum treatment of the dust col¬ 
lapse applicable to ultra-light particles in Schwarzschild 
coordinates. The conclusions follow in Section V. 


I. THE OPPENHEIMER-SNYDER MODEL 


In General Relativity, a first approximation to the 
exterior space-time of any star, planet or black hole 
is a spherically symmetric space-time modelled by the 
Schwarzschild metric. This is a consequence of Birkhoff’s 
theorem [3]. The Schwarzschild line element thus takes 
the usual form 


ds 2 


where 


2 M\ o dr 2 

- ) dt 2 -\ -- 

r 1 - 2 Mr 




( 1 ) 


dVt 2 = dd 2 + sin 2 9d(f) 2 . (2) 


Throughout this paper we use geometric units with G = 

c = 1 . 

The Oppenheimer-Snyder (O-S) model follows the col¬ 
lapse of a star that is idealized as a dust sphere with 


uniform density and zero pressure from the perspective 
of an observer located on the surface of the star. The 
motion of the collapsing surface initially at radius r i can 
be parametrized by 


= 7T (1 + cos V) 




t = 2 Mlog 


\Jrt/2M — 1 — tan 77/2 


\Jtij2M — 1 + tan 77/2 




v ■ 

v + ^(77 +sin 77 ) 


( 3 ) 

( 4 ) 

( 5 ) 


These equations correspond to Eqs. (32.10a) to (32.10c) 
of [3]. The star collapses to a singularity in finite proper 
time. However, it takes an infinite Schwarzschild t to 
reach the apparent horizon at r = 2 M and thus an ex¬ 
ternal observer will never see the star passing its gravi¬ 
tational radius (r = 2 M). 

In [3j the stellar surface is considered to be initially 
at rest. Here we integrate the equations of motion for 
configurations with all possible initial velocities. 


II. CLASSICAL ACTION AND RADIAL 
GEODESICS IN SCHWARZSCHILD 
COORDINATES 


The relativistic action S and the Lagrangian C for this 
system are 


S c 1 = —771 I dr = —m I dt\! ( 1 — -- ) — 


J 1 - 2 M/r ’ 

( 6 ) 

where r is the proper time, and r = dr/dt. The La¬ 
grangian is 


C = -mi / 1 


2 M\ 


r J 1 — 2 M/r 


The momentum 

dC 

P= ^f = 


mr 


'(1 - 2 M/r) |^(1 - 2 M/r) - f 2 
and the Hamiltonian is 


( 7 ) 


( 8 ) 


2 M 


H = pr — C = 1- \lm 2 +p 2 [ 1-— J, (9) 


2 M 




which reduces to the free particle Hamiltonian in the 
Af —^ 0 limit. The equation of motion derived from the 
Euler-Lagrange equations is 


3 M / dr \ 
r 2 \ dt J 



d 2 r M / 2 M 

dt 2 + r 2 \ r 


= 0 , 
( 10 ) 
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which can be integrated to 

2 




2 M . 

Ci H-(1 — Ci) 

r 


From this the classical paths are determined 


[tf r r f 

/ dt = tf—ti = ± dr- 
J t; J Tj v 


where 


x r {r — 2 M) 
= \J c±r 2 + 2 Mr(l — ci). 


( 11 ) 

( 12 ) 

(13) 


The + sign represents motions from ?y to 77 > ?y, and 
— sign represents the star collapsing from 7 to 77 < tv 
T he parameter ci can be used to separate the space-time 
regions 


ds 2 > 0 
ds 2 = 0 
ds 2 < 0 


ci > 1 outside the light cone (14) 
Ci = 1 on the light cone 
ci < 1 inside the light cone 

(15) 


In the free particle case (M —► 0 in Eq. (Ill), ci = r 2 
represents the velocity squared of the particle. 


Stellar Surface at Rest Limit 


When the surface of the star is initially at rest, initial 
velocity 


= 0 . 


From Eqs. this corresponds to 

2 M 

' r, - 2 M ' 


ci = 


(16) 


We can recover our equations of motion for this ci value 
from the O-S model. We proceed by taking the derivative 
of Eq. ([3| and Eq. Q w.r.t. 77 

dr r.j 


drj 

dt 

drj 


sm 77 


= ~\/r{ri - r) 


(17) 


.mvIT 

cos 4 (?7/2) 

V 2 M 

_( ri/2M ) cos 2 (t7/2) — 1 

2 

- 1 , 

2 M V 2 M 


where we used that r = n cos 2 rj/2. We then divide dr / dr] 
by dt/dr], and use r* = 2M{c\ — l)/ci from Eq. (16) to 
recover the equation of motion 


%=-- — ^—V c i r2 + 2 Mr(l - ci). (18) 

dt r z 

The negative sign — fits in with our convention for a 
collapsing star. While some paths will have low proba¬ 
bility, a path integral approach that includes quantum 
effects requires the consideration of space-like, time-like 
and light-like paths with all possible initial velocities. 


A. Classical Paths in Schwarzschild Coordinates 


Schwarzschild coordinates allow us to explore the 
space-time outside r = 2 M. The possible paths taken 
by a particle on the surface of a collapsing star are writ¬ 
ten explicitly in Table [I] Every point {rj,tf) may be 
reached either via a direct path or via an indirect path. 
Table |H] contains the boundary regions that determine 
whether a path is direct or indirect, and whether it lies 
inside or outside the light cone. 

The table uses 


x(r s ) = x s = \Jc x r 2 + 2Mr s (l - Ci), 


(19) 


0s = yfcx, + Cir s + M( 1 - Cl), 


( 20 ) 


U s = 


Cl 


2 M log(r s - 2 M) - 2 M log 


{r s +x s f 


( 21 ) 

where s = i, /, a. 

The direct paths include a logarithmic term when 
Ci > 0 and an arc sin term when Ci < 0. The indirect 
paths occur for r t in the 63 — 64 region. They each contain 
a zero velocity point where the path turns around. Math¬ 
ematically, the velocity passes through zero ( dr/dt = 0 
at some radius r a ) when x(r a ) = 0, which corresponds to 
ci = —2 M/(r a — 2 M). The indirect paths reach the final 
point ( r/,tf) after going through a zero velocity point at 
r a (t a < tf) with r a > rf and r a > rt. The paths from r, 
to r a correspond to the positive sign in Eq. (12), while 


the path from r a to rf corresponds to the negative sign. 

To understand these paths we make an analogy with 
the vertical motion of a ball under gravity (See Fig. [I]). 
This analogy arises naturally due to the correspondence 
between the Friedmann equation and the equation of mo¬ 
tion for a test particle in the gravitational field of a point 
mass in Newtonian gravity. It can be made only in the 
Schwarzschild case because then the interior of the star 
can be represented by an FRW model, which provides 
the link to Newtonian gravity. 

Consider a ball traveling from point (initial radius) 
to point rf (final radius) with r/ below r^. It can either 
go down directly from r* to rf (direct path) or it can go 
up from r'i , reach its highest point, where its velocity will 
be zero, and then fall back down to rj (indirect path). 
This is shown schematically in Fig. [lja). Analogously, if 
a ball has to travel from r l to r/ with rf above r t . it can 
go directly from r* to rf or it can pass above r/, reach 
its zero velocity point above rf and return down to rf. 
This case can be seen in Fig. So- The middle figure 
(Fig. [ljb)) shows the special case when n = 77. Then 
the ball can either stay put or it can be thrown up, reach 
its higher point, and fall back down. Each of these paths 
takes a different amount of time to complete. For a fixed 
travel time At, the two points r, and rf are connected by 
a single path that is either a direct path or a turn-around 
path. At the turning point of a turn-around path (also 
its highest point) the speed of the particle is always zero. 
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Range of 14 ci 


Equations of Motion 


t f -U=a f -ai + M^ ^log(^) 

t f -t i =a f -a i + /fgA- 1 [arcsin “ *) “ arcsin ( m~(i-c1)T “ *)] 

tf-U= 2a a -aif-an + M [2 arcsin ( 1 ) - arcsin (jg/g - l) - arcsin - l)] 

ti- tf = a f -ai+ M c f^ [arcsin ( g/g, - l) - arcsin ( y/l’g - l)] 
ti- tj =a / -Q i + M^^log(|f-) _ 


2 M - b 2 
b 2 - 63 

63 - 64 

64 - 65 

65 — 00 


ci > 0 


2 .M n 

r f -2M < C 1 < 0 

2 M 

r a —2M 

2M / _ / n 
64 - 2 AI < Cl < 0 

ci > 0 


TABLE I: Closed-form solution to Eq. (12) in Schwarzschild coordinates. At each region boundary bj (j = 2 — 5) the 
equations converge to the form given in Table 0 All discontinuities cancel resulting in smooth transitions between 

regions without singularities. 


Cl 


Regional Boundaries 


tf~U=r f 


V2 ( 3/2 

■ (/ 


2A1 

' 64-2M 


^ ^ “ 3CM 

2Mlog(l-^) - ^+2Mlog 

2 M j 


tf ti 


t,-t / =^- 2 Mlog(l-f 


\ Ol —ZiW / 

- & 3/2 ) + 2\/2M (Jrj - /fe) + 2Mlog 

V) - TT + 2M1 °S + M fil 

(3ci-1) 


( i + i/vf 


f )(‘ 


jrMl [ arcsin ( X ) “ arcsin ( «M)' 


t — tf — ^ ( r 3 ' 

~ 3 Vm V / 

U-tf = r f -b 6 + 2M log 


ry y ci 1 & 63(63 — 2 M) civ/—cl L ' ' 

(1 - 2 M) - 2 Mlog [ r % + :^) ] + [- arcsin ( 1 ) + 

/Si j C + /23T' 

6 3/2 ) + 2\[2M (y/rj - s/b$ + 2M log ) V gg -— 

fin <r( r f~ 2M \ 


( M(l-g) X )] 


arcsin ( 1 ) + arcsin ( yg/gy - l) 

/ O A/f \ 


TABLE II: Region boundaries for the equation of motion (Eq. ©)■ The first column is the border value of ry for 
each region given the r f and ci values. It is calculated using the equation of motion in the third column. The light 

cone boundary corresponds to ci = 1 . 


If we keep the time interval constant, we find that for a 
given Tf there is a point r* = 64 > r/, where the velocity 
dr/dt at 64 = 0. Similarly, there is a point = 63 < r/ 
such that the velocity dr/dt at rf = 0. For any point ?y 
between b 3 and 64 , the only way a particle can reach rf 
in the same time At is if it goes through a turning point 
r a > rf and r a > n. The particle cannot reach rf>Ti 
if it’s velocity vanishes at any point between r,; and rf. 

The motion towards or away from r = 2 M can be com¬ 
pared to the vertical motion of a ball under gravity when 
thrown straight down or up, respectively. Classically, 
the star will always collapse to a black hole. Quantum 
mechanically, we include all possible initial velocities of 
the surface of the star towards and away from 2 M. The 
wavefunction of a particle on the surface of the star de¬ 
termines the spread of the trajectory of the surface of the 
star itself. 

Escape Velocity. For a particle of mass m that is 
thrown up from the surface of the Earth, an initial up¬ 
ward velocity Escape > 11.2 km/sec ensures that the par¬ 
ticle escapes Earth gravity. At this speed 

GM E m _ mvl 

Re 2 ’ ( j 

where Me, and R E are the Earth mass and radius. 


We continue with the ball analogy to estimate the es¬ 
cape velocity for our particle. In our case, R E is replaced 
by the initial radius of the star rj, and Vi = dr/dr. The 
kinetic energy of the particle on the surface of the star 
equals its potential energy when 


Mm 

n 


m 

2 ~ 


Here 




dr 

dr 

2 M 
r 


-1 


dr 

dt 


with dr/dt = r given by Eq. @. This results in 


Mm me i Mm, 

-= +- (1 — ci), 

Tj 2 Vn 


(23) 


(24) 


(25) 


whose solution is Ci = 0. So, while all outgoing paths 
with ci < 0 will turn-around, outgoing paths with ci > 0 
will not have turning points. When ci = 0, the turn¬ 
ing point (dr/dr\ Ta = 0 , x a = 0 ) occurs at r a —> 00 as 
expected for an escaped particle. 

The different paths (Table ^ and their various re¬ 
gional boundaries (Table [Tl| are summarized below for 
a given (77, t/) outside r = 2 M: 
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FIG. 1: Schematic of paths, 
(a) ry < n 


• Ci = —2M/(rf — 2M): ry = 63 < r/. The mini¬ 
mum value of ?y where ry is reached directly with 
aj(ry) = 0 and null final velocity 


turn around radius 



(b) 77 = n 

turn around radius 



(c) 77 > ry 



initial radius 


• Ci = 1 : the paths arc light-like ( ds 2 = 0 ) with ry = 
b\ < ry and ry = be > r / . The paths with Ci > 1 
are space-like, while those with Ci < 1 are time¬ 
like. 

• Ci = 0: Ti = b ‘2 < ry and ry = 65 > 77 . This 
path defines the escape velocity. All paths with 
Ci > 0 that are outgoing will not reverse, allowing 
the particle to escape the gravity of the star. 


Paths with ry < b 3 are direct. Paths with ry be¬ 
tween 63 and 64 are indirect. 

• Ci = -2M/(b 4 — 2M): ry = 64 > 77 . This path 
starts with zero initial velocity. All paths with ry > 
64 are direct. 


III. CLASSICAL PATHS IN KRUSKAL 
COORDINATES 


Our study is the first to consider space-like paths in 
Kruskal coordinates for a particle on the surface of a col¬ 
lapsing star. They have not been explicitly considered 
before because they are believed to be of low probability. 
We conjecture that they play an important role in quan¬ 
tum mechanical collapse in the same way the probability 
of passing through the potential barier is important in 
tunneling. Classically, it can never happen, and yet this 
probability cannot be ignored in quantum mechanics. 

In this section we derive and discuss the closed-form 
solutions to the equation of motion for the space-like, 
light-like and time-like paths in Kruskal coordinates. We 
find that space-like geodesics have the interesting prop¬ 
erty that they can turn around outside r = 2 M and 
move back in time. Unlike the Schwarzschilcl t and r, 
which are time and space coordinates respectively out¬ 
side r = 2 M but switch roles for r < 2 M, the Kruskal 
coordinate v is always a time coordinate and the coordi¬ 
nate u is always a space coordinate. The relation between 
the Schwarzschild r and t and the Kruskal variables u and 
v is given by 


u = 


with 


and 


, / JL _ 1 p r/4M CQS u _±_ 

V 2 M iC CObI1 4 M 

r > 2 M 

, /l_ r_ r/4M • k _±_ 

V 1 2M e blIH1 4 M 

r < 2 M 

. / JL _i p r/4 M s ; n u J_ 

V 2 M 1C blliIi 4 M 

r > 2M 

Vi - m er/4M cosh tm 

r < 2 M 

u 2 - v 2 = (— - 1 ) e 1 

72 M 

\2M ) 



v , t „,, 

— = tanh —— : , r > 2 M 

u 4 M 

u , t „,, 

— = tanh-, r < 2 M. 

v 4 M 


(26) 

(27) 


Thus the line element in Kruskal coordinates takes the 
form 

32 M 3 


ds" = 


-e~ r / 2M (du 2 -dv 2 )+r 2 dn 2 . 


(28) 
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The singularity at r = 0 occurs at u 2 — v 2 = — 1. In 
this paper we consider only the first quadrant u > 0 and 
v > 0. The horizon r = 2M occurs at u = v. Outside 
the horizon u > v, whereas inside the horizon u < v. 

The classical equation of motion 


u 

1 — u 2 


= 1 - 


4M- 


can be derived from the Euler-Lagrange equations 


d_ DC _ dC 
dv du du 


(29) 


(30) 


for the Lagrangian 


C(u, u, v ) = -4 MJ (1 U P 2 {1 * M/r) - (31) 

y [u z — v z ) 

Note that here u = du/dv. 

The equation of motion can be reduced to 


du v/u±x r /r 
dv 1 ± ( v/u)(x r /r ) ’ 


(32) 


where x r is given by Eq. Direct paths from ry 

to r/ > Ti have the “+”-sign in the numerator and de¬ 
nominator with du/dv > 0. Direct paths from r, to rj 
with Tf < r,j will have the “—’’-sign in the numerator and 
denominator. 

The classical paths describing the evolution of the star 
outside r = 2 M have been discussed in Sec. EEa] in 
Schwarzschild coordinates. It was shown that there were 
direct space-like and light-like paths as well as direct and 
indirect time-like paths. Any two points in (r, t) were 
connected by a unique classical path. 

What is new in Kruskal space-time are the turning 
points in the u — v plane for all space-like, light-like and 
time-like paths. Spacelike paths (ci > 1, \du/dv\ > 1) 
turn in time at points when dv/du = 0 , whereas time¬ 
like paths (ci < 1 , \du/dv\ < 1 ) turn in space when 
du/dv = 0. Additionally, there are indirect space-like 
and light-like paths that have turning points in r inside 
the apparent horizon (here we refer to paths as being in¬ 
direct when they turn in r; the light-like paths turn at 
r = 0.) In contrast to the time-like paths, the in-going 
space-like paths are not unique. A point outside r = 2AI 
might be connected to a point inside the apparent hori¬ 
zon via a direct and an indirect space-like path or via 
two indirect space-like paths. 

We determine the regions (Table |III[ ) and the light- 
cone boundaries (Table Rvf that delimit the different 
kinds of paths that start at ( Ui , Vi = 0 ) outside the hori¬ 
zon and reach (uf,Vf) inside the apparent horizon. Out¬ 
side the horizon, the light cone boundary values of ‘zq’ 
for a given Uf and vr can be determined by replacing 
tf — ti in Table [n] by tanh _ 1 (u//M/). The bound¬ 
ary values of n (Vi = 0 ) with 77 < 2 M, are calculated 


from the equations in Table IV Note that Uf is a func¬ 
tion of Vf and 77 and Ui depends on zy, which is typ¬ 
ically taken to be zero, and Ty. Thus if one finds r, 
this determines zq. The point r,; = b' 4 corresponds to 
Xi = 0 with ci = —2M/{b\ — 2M). The Ci = 0 and 
ci = 1 (lightlike) paths that reach points (uf,Vf) in¬ 
side the horizon are used to determine the r, values b' 5 
and b' 56 respectively. The 7 y = b' 6 value corresponds to 
Xf = 0 with Ci = — 2 M/(rf — 2M) (turning point at 
r = rf < 2M, space-like path, ci > 1). The value of 
n = b' 7 is to be determined by using ci = - 2 M/(r a - 2 M) 
(x a = 0,r a < r/ < 2 M) in the equation of motion. The 
two unknowns b 7 and r a are determined by using the 
equation of motion and its partial derivative with respect 
to Ci. Both equations are given in Table |IV| In the re¬ 
gion outside b' 7 (> b' 7 ), there exist no classical paths to 
(u f ,Vf). 

In Fig. [2ja) and (b) we trace direct time-like, di¬ 
rect space-like, and indirect space-like trajectories that 
start at z;j = 0 and pass through a fixed point P(zzp = 
9.995, z;p = 10.0) with rp ss 0.8 M, which is located in¬ 
side the apparent horizon. Path 1 and path 2 are time- 
like, and originate at different points outside r = 2 M. 
Path 1 (Ui = 2.59, ci = —1.05) lies in the b\ — 6' 5 re¬ 
gion of Table III Path 2 («., = 16.15, cl = 0.8) lies in 


the 65 — b' 56 region, which contains time-like paths with 
a turning point in the u — v plane ( du/dv = 0 ) that lies 
outside the apparent horizon. 

The next region is b' 56 — 6 g, where each point out¬ 
side the horizon is connected to points inside the hori¬ 
zon by two space-like paths: one direct and one indi¬ 
rect. Path 3a and path 3b originate at the same 
Ui = 27.6 (r* = 10.5 M). Path 3a reaches P directly. 
However, after having passed through P, it turns at 
r a = 0.48 M before reaching the horizon like all space- 
like paths. Path 3b is indirect turning at r a = 0.3 M 
( u a = 11.3, v a = 11.35), and then reaching P on its way 
back towards the apparent horizon. Both paths also have 
a v — u turning point outside the apparent horizon where 
dv/du = 0. After having passed through its v — u turning 
point, each path moves back in time (the space coordi¬ 
nate u continuously decreases, while the time coordinate 
v increases up to the turning point and then decreases). 

In the r,; = b' 6 — b 7 region there are two indirect 
paths that connect the same point outside the horizon 
to a point inside r = 2 M. Paths 4a and 4b con¬ 
nect the point of zq = 46.3 to P. They each have a 
dv/du = 0 turning point outside the horizon, and also 
turn in r inside the horizon before reaching P. Path 
4a turns at r a = 0.77147 < rp and path 4b turns at 
r a = 0.79969 < rp. It can be seen that the paths are 
very close together. As Uj increases, the paths in this 
region become progressively closer until they merge at 
rj = b 7 . Path 5 shows the single indirect path that orig¬ 
inates at Ui = 46.7 (r,; = b 7 = 12.13 M). Beyond this 
point, there are no real solutions and hence no way to 
reach point P. 

Fig. § c) displays the c\ values of space-like paths 
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Range of ri ci 


Equations of Motion for Paths with r, > 2 M and 77 < 2M. 


-AM tanh 



ctf — an + |^arcsin ^ 

a f -ai + M^A log (ft) 

= 2« a - a/-ai + M^^Iog(^ 


l) — arcsin ( 


64 - &5 

6 5 — &6 
&56 — K 


2A4 ^ n 

b' 4 -2M < C 1 < 0 


l r / 


o 


< Cl < - 
Cl = 


2 M 


^ 4M t anh- 


M(l-ci) 


Af(l-ci) 


-2M 


-4A7 tanh 


L) 

<f) 


TABLE III: The equation of motion with rj < 2 M. At the light cone boundaries they converge to the form given in 
Table IV All discontinuities cancel each other and there are no singularities anywhere in the equations of motion. 



TABLE IV: Equations of motion at light cone boundaries for paths with 77 < 2 M. Outside 2 M the equations for 
Table Imold with tf — ti replaced by AM tanh” 1 (17/17). Like before, the value ci = 1 corresponds to a light-like 
path. Note that b 7 is the caustic point where both the path /(— 2 M/(r a — 2 M), u,v) and its partial derivative 

df/dc 1 vanish. 


reaching point P (up = 9.995, up = 10.0) from m > 46.1 
at Vi = 0. At first, for each 14 in the figure there are two 
ci values with which the final point P can be reached. 
The figure clearly shows the roots (ci values) coming 
closer and closer together until they merge at the caus¬ 
tic point r = b 7 . Beyond this point there are no paths 
that reach point P. Extensions to this work to determine 
the Kruskal-WKB wavefunction will require an analysis 
of the caustic at r = b 7 to remove any potential diver¬ 
gences. 

The indirect paths connecting {u l . 7 ) with 7 > 2 M 
to ( up,vp ) with rp < 2 M penetrate the horizon deeper 
than rp turning at r a < rp before reaching rp (e.g., see 
the green path in Fig. [2](a)). After turning, all paths 
must continue towards the apparent horizon reaching it 
at u = v = 0. These paths move back in time from a 
point outside the horizon, where dv/du = 0 , after hav¬ 
ing traveled to that point forward in Kruskal time from 
= 0 ) with 7 > 2 M. 

As discussed before, the turning points in r corre¬ 
spond to x r = 0 (in the Kruskal paths x r appears in 
Eq. (58)). They occur at r = r a when x a = 0, where 
cl = —2M/(r a — 2 M). Since Ci > 1 for space-like paths, 
clearly r a < 2 M (the turning points are inside the hori¬ 
zon). For any r on a path with r a as the turning point, 
x r can be written as yj2Mr(r — r a )/(2M — r a ). To en¬ 
sure the non-negativity of the term under the square root 


all points on this path must have r > r a . By the same 
token, for time-like paths r a > 2 M and r < r Q . This case 
was described in the Schwarzschild analysis. 

All space-like paths that reach points inside the horizon 
subsequently head towards the horizon r = 2 M (u = v). 
Classically, they end at Uf = 77 = 0 . However, quantum 
mechanically there may be paths close to the classical 
path that bring information from within the black hole 
horizon to the outside. 


IV. QUANTUM TREATMENT IN 
SCHWARZSCHILD COORDINATES 

We consider a collapsing cloud of ultra-light particles 
with rriM ~ h. Such particles are not localized and thus 
neither is the surface of the star. The position of a par¬ 
ticle on the surface of the star is then approximated by 
a wavefunction. In a relativistic path integral approach, 
this wavefunction can be computed in the WKB approx¬ 
imation via an integral over the classical action | 23 ] • It is 
necessary to include configurations will all possible initial 
velocities. For each point ( rf,tf ) all classical paths de¬ 
rived above are needed. In this paper, we only perform 
the WKB analysis in Schwarzschild coordinates, where 
we can only use the paths to 77 > 2 M. 













































TABLE V: The classical action in various regions. The + sign is chosen when rj > r*, and the — sign when r/ < r,;. 


A. WKB Approximation: Closed-form Propagator, 
Numerical Wavefunction 


In order to study the wavefunction we use a WKB 
approximation of the propagator 


G(r, t; Ti, U) = [ Vrexp[iS/h\, (33) 

Jc 

where S is the action associated with each path. The 
set of paths C include space-like, time-like and light-like 
paths. 

The WKB approximation involves expanding about 
the classical action 5' c i to include paths that slightly de¬ 
viate from the classical paths. The propagator is dom¬ 
inated by paths near the classical trajectory between 
(■ ri,U ) and (r,f) and is approximated by the WKB ex¬ 
pression [24] 

Gwkb (r, t; n, U) = Aft exp (i : X ) . (34) 


The WKB wave function is obtained from the integration 
of the propagator 

pOO 

^wkb i r f,tf)= / GwKB(’rf,tf,r i ,ti)'S'(r i ,ti)dr i . 

Jo 

(35) 

The initial wavefunction that describes a particle on 
the surface of a star of radius r, is taken to be a Gaussian 
centered about r c 


n,ti = 0) = exp 



(36) 


with r c » 2 M. 

The action from Eq. © is rewritten using Eqs. 
and (13) as 


r( r t^f) r 

S c i(ri,tf,rf,tf) =±imy/c 1 -l / dr ^T\ ( 37 ) 

J(ri ,ti) X[T) 

for direct paths connecting (r,,tj) to (r/,f/). The + sign 
corresponds to rj > and the — sign corresponds to 

Tf < Ti. 

Indirect paths connecting (n, ti) to (r/, tf) through the 
turning point ( r a ,t a ) with r a > r.j and r a > rj are de¬ 
scribed by the classical action 


S d = im\J C\ — 1 


x(r) 


L 


(Vf,tf) r 

dr—- 

(r a ,t a ) x \ r ) 


(38) 

In Schwarzschild coordinates, the wavefunction is con¬ 
tinuous across all regions outside r = 2 M. The presence 
of the r = 2 M coordinate singularity prohibits the inclu¬ 
sion of paths at or beyond the horizon region. Closed- 
form expressions for S c \ and d 2 S c \/(dr.idrf) are given in 
Table \V\ Both functions are continuous across all re¬ 
gions outside r = 2 M. Once ci is known for every point 
(r, t), we have all the ingredients to find the propagator 
in Eq. (34) using the S c \, and d 2 S c \/dridrf expressions 
from Table [V] The propagator is then integrated to find 
the wavefunction. 
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the wavefunction in the WKB approximation converges 
to the relativistic Schrodinger solution, which is the ex¬ 
act solution for this case. In [7] some disagreement was 
observed, but while they attempted a physical interpre¬ 
tation, it is likely due to numerical error. 

When M —> 0, the Hamiltonian from ([9]) reduces to 

H = a/p 2 + to 2 , (39) 


where p is now the momentum operator. The Hamilto¬ 
nian is nonpolynomial in p. The square root term thus 
corresponds to a non-local momentum operator that is 
interpreted as acting on any wavefunction [7] 

/ OO 

dke ikx ct)(k,t). (40) 

-OO 


to give 

H^(x,t) 



dke zkx \/h 2 k 2 + m 2 cj)(k, t). 


(41) 


This integrates to 


4>{k, t) = A(k, xq) exp 




(42) 


When At = 0, we recover the original wavefunction, 
which is chosen to be a Gaussian centred around the ori¬ 
gin, 


^(x 0 ,t 0 ) = A/" exp 


m 2 Xg \ 

h 2 ) 


and so 


4'CM) 


W 

27r 


/ oo /»oo 

dk / dxo exp (— m 2 xl/h 2 ) (43) 

-oo J — oo 


X 


exp (ikAx) exp 




FIG. 2: (a) v vs r and (b) v vs u are shown for 
space-like and time-like paths with Vi = 0 that pass 
through P(up = 9.995, up = 10.0) of rp = 0.8, a point 
located inside the apparent horizon, (c) The C\ values 
for indirect space-like paths in the region between U 6 
and b’-j from which point P can be reached are displayed 
for each Ui. No paths exist beyond r = b' 7 , which is 
where the two roots merge into one. 


B. Relativistic Schrodinger Solution 

1. Free Particle 

The M —> 0 limit describes a relativistic free particle 
that is no longer confined by the gravity of the sphere 
of dust. This case was first studied by Redmount and 
Suen in [7.. We revisit this limit to investigate whether 


The normalization \Af\ 2 = m\/2/(fty/ir) is time indepen¬ 
dent. 

In terms of the propagator 


/ OO 

dz 0 G(z/,f/;x 0 ,to)^(2:o,io), (44) 

-OO 


where 


/ °° ( ik 

— exp (ikAx) 

-OO 


x exp 


— —~V h 2 k 2 + m 2 

h 


= ! rdAl±£) Kl ( n ayy B ).(45) 

e->o+ 7m v A e 


Here K\ is the modified Bessel function and X e = 
Ax 2 + (iAt + e) 2 . When At = 0 the propagator reduces 
to S(Ax) as expected. 
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FIG. 3: Re[^], Im[’F] and |'F| 2 are shown for the free 
particle case at t = lOh/m as a function of Xf, which 
also has units of h/m. The solid lines show the 
Schrodinger solution, while dotted lines and dashed lines 
show the WKB approximation for e = 0.05, e = 0.005 
respectively. It can be seen that as e —► 0 the WKB 
approximation converges to the Schrodinger solution. 


The WKB propagator for the relativistic free particle 
is then computed using Eq. (34) >7j 


Gwkb = 


m (iAt + e) 


27r 


n a; 


3/2 


exp 


m\ l J 2 

h 


(46) 


where WKB wavefunction is given by 




WKB 


{X,t) = [ 

J —c 


dx 0 G W KB(x, t; Xo, t 0 )T(x 0 , t 0 ). 

i 

(47) 

are evaluated 


The integrals from Eq. (44) and Eq. ( |4 
numerically. The results are shown in Fig. [3] It can be 
seen that the wavefunction in the WKB approximation 
converges to the Schrodinger solution in all parts of the 
light cone. When no star M is present, the relativistic 
Schrodinger is exact as originally discussed in Redmount 
and Suen [7j. However, when they performed the same 
numerical comparison, they found some disagreement be¬ 
tween the two solutions. This we believe to be due to 
numerical error of [7]. We observe convergence (see Fig. 
i) of the wavefunction in the WKB approximation to 
the Schrodinger solution in the M —> 0 limit. 


FIG. 4: The time evolution of WKB |'3/(rj,fy)| 2 is 
shown from t = 5 M to t = 45 M in the case of 
mM/h = 1. We can see that the star collapses to a 
black hole. 


where p is the momentum operator and 

B = B(M,r) = 1 — -—. (49) 

r 

The whole wavefunction can then be written as 

T(r,t) = J driG(r,t-,ri,ti)^(ri,Ti), (50) 

where the propagator is given by 

/ °° dk 

— exp(zfcAr) (51) 

-oo 2n 

-i -— \J Ti 2 k 2 B 2 + m 2 B 

n 


x exp 


Like before, the propagator can then be integrated ex¬ 
actly to obtain 

G(x, t; Xi, ti) = lim m (^ + e ) (mB 1 / 2 \ l J 2 /%) (52) 

nhB l G\l /2 

where K\ is the modified Bessel function and 


K = 


A r 

13 


+ (i At + e) 2 . 


(53) 


2. Non-zero mass 

While in the free particle case the Schrodinger solution 
was exact, when M ^ 0 it becomes a rough approxima¬ 
tion, which fails as we approach r = 2 M. Comparing the 
two solutions is still instructive. If we follow the same 
procedure as for the free particle above, and re-write the 
Hamiltonian from Eq. as 

H = yj B 2 p 2 + Bm 2 , (48) 


When r 2 M the propagator vanishes. Thus this so¬ 
lution is inaccurate at late times and cannot model the 
final stages of the collapse of the dust sphere. 


C. Numerical Results 


We integrate Eq. (35) numerically for the Gaussian 
wavefunction from Eq. (36) centered about r c = 30 M. 


Classically, the star collapses to r = 2 M in infinite 
Schwarzschild time. However, quantum mechanically 
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FIG. 5: (a) I'E'fV/, t /)| 2 at t = 5 M as a function of r/ 
for mM/h = 4.0, mM/h = 1, mM/Ti = 0.5 and 
mM/h = 0.25. It can be seen that that lower masses 
behave more quantum mechanically, and are less likely 
to collapse to a black hole, (b) |'F(r/, f /)| 2 at t = 10 M 
as a function of rf for the same masses. 



FIG. 6 : mM/h = 1 Schrodinger (solid line) and WKB 
comparison (dotted line). 


there is a possibility for the star expanding and dispers¬ 
ing as well. Fig. [4] shows the time evolution from t = 0 
to t = 45 M of \'$’wKB( r fitf)\ 2 as a function of rf. An 
ingoing peak and an outgoing peak appear with the am¬ 
plitude of the ingoing peak always remaining larger than 
that of the outgoing peak for the course of the evolution, 
which is expected for a collapsing star. For numerical 
computations we scale all variables to be dimensionless. 


The only parameter is mM/h. We note that Fig. [4] is 
qualitatively similar to the WKB results of pL2l . 

Fig. [5](a) and Fig. [5](b) show the ingoing and outgo¬ 
ing peaks of \^wKB(f, tf )\ 2 for different values of mM/h 
at times of tf = 5 M and tf = 10M, respectively. For a 
given particle number N = M/m , the speed of the col¬ 
lapse increases slightly with the particle mass m as does 
the value of the ingoing peak. This is consistent with 
our expectation that for higher particle mass the star is 
more classical. For lower masses, the star resists collapse 
longer until probability of dispersion exceeds the proba¬ 
bility of collapse. 

In Fig. [ 6 ] we compare the WKB and Schrodinger so¬ 
lutions. It can be seen that the radial position of the 
ingoing peak drifts out of step at late times, while the 
outgoing peak continues to evolve at about the same po¬ 
sition. As the star approaches rf = 2 M, the ingoing 
peak amplitude for the WKB approximation increases 
developing a numerical singularity that indicates the for¬ 
mation of an apparent horizon, while the ingoing peak 
amplitude of the Schrodinger solution decreases since the 
propagator in Eq. (52) vanishes at r/ —> 2 M. The lat¬ 
ter behaviour reinforces that the relativistic Schrodinger 
equation is not the correct representation for a particle 
in the gravitational field of a collapsing star. 


V. CONCLUSION 

Closed form solutions for the classical paths taken by 
a particle on the surface of a collapsing star have been 
determined for all initial configurations in Schwarzschild 
and Kruskal coordinates. 

We found that 

(i) all time-like paths are unique. For a given time in¬ 
terval, a path between an initial and final point can 
be either direct or indirect, where it turns around 
in space. Thus some particles that initially move 
away from the star, can return and contribute to 
the collapse. 

(ii) all space-like paths are unique outside r = 2 M. 

(iii) Kruskal space-like paths can turn around in 
Kruskal time, but cannot turn in Kruskal space. 
Multiple paths can connect a given point outside 
the horizon to a final point that lies inside r = 2 M. 

(iv) Classical Kruskal space-like paths connect points 
outside r = 2 M with points inside r = 2 M up to 
a critical value of r = r a . Spacelike paths from 
(r i , Vi = 0) to any given point r inside the horizon, 
where r a < r < 2M, are non-unique (two paths 
exist). These paths get less and less separated and 
merge into a single space-like path at r = r a and 
any point r < r a is unreachable from n. Upon 
reaching the critical value r = r a , the space-like 
paths turn back towards r = 2M , reaching it at u = 
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v = 0. Therefore, classically, no information from 
r < 2 M can exit a black hole. However, by taking 
into account paths close to the classical paths, one 
might be able to extract information from r < 2 M. 

If we make an analogy to tunneling, the particle 
has a low probability of passing through a poten¬ 
tial barier. Classically, it never happens, and yet 
this probability cannot be ignored in quantum me¬ 
chanics. Similarly, it may be that the low proba¬ 
bility space-like paths will play a crucial role in our 
understanding of quatum mechanical collapse. 

The collapse of a self-gravitating star was next modeled 
as a ball of dust using a WKB approximation. We inte¬ 
grated around the classical paths with all possible initial 
velocities to include quantum effects. This extended the 
analysis of [7] from the relativistic free particle to the case 
of non-trivial gravity. In practice, quantum mechanical 
effects are important in macroscopic stellar collapse when 
the sphere is composed of ultra-light particles. A number 
of ultralight dark matter particles have been proposed, 
which could physically motivate such stars. 

The evolution of the wavefunction of a particle on 
the surface of a collapsing star was followed numeri¬ 
cally. We showed that in the case of a star collapsing 
with zero initial velocity, our path equations reduce to 
the Oppenheimer-Schneider equations of motion. In the 
M —> 0 limit, the relativistic free particle wavefunction of 
[7] is obtained, and the WKB and relativistic Schrodinger 
solutions match. Since in this limit the Schrodinger so¬ 
lution is exact [253, the covergence of the WKB to the 
Schrodinger solution enforces the validity of the WKB 
approximation. 

Our results for this self-gravitating collapse are sum¬ 
marized as follows: 

(i) The wavefunction representing a particle on the 
surface of a collapsing star typically exhibits an 
outgoing and an ingoing component, where the for¬ 
mer contains the probability that the star will dis¬ 
perse and the latter the collapse probability. For a 
given particle number N = M/m, we find that the 
rate at which collapse occurs increases with particle 
mass. This is consistent with the expectation that 
for higher particle mass the star is more classical. 
As the particle mass is lowered, the star resists col¬ 
lapse until the probability that it disperses exceeds 
the collapse probability. Note that some part of the 
star always disperses even when a black hole forms. 

(ii) In the case of the collapsing star the relativistic 
Schrodinger solution is not a good approximation 
for the wave function. On comparing the WKB 
and relativistic Schrodinger solutions, we find that 
the ingoing wavefunction gets more and more out 
of step at late times. The outgoing component of 
the wavefunction for the two solutions is in bet¬ 
ter agreement because it is further from the coor¬ 
dinate singularity at r = 2 M. As expected, the 


Schrodinger and WKB solutions are out of step at 
early times when WKB solution is more inaccurate, 
they come closer together at intermediate times, 
and fall out of step again at late times, when the 
Schrodinger approximation fails to model the sin¬ 
gularity formation. The presence of the coordinate 
singularity at r = 2 M, motivate a potential inves¬ 
tigation beyond r = 2 M via the WKB approxi¬ 
mation deployed in Kruskal coordinates where no 
analytical relativistic Schrodinger solution is avail¬ 
able. 

The traditional Oppenheimer Schneider model de¬ 
scribes the zero pressure collapse of a star starting from 
rest. We use a path integral approach that requires in¬ 
cluding paths of all velocities. For a particle on the sur¬ 
face of a star, we explicitly parametrically describe all 
such paths of different initial velocities. This is done in 
Schwarzschild coordinates until r = 2 M and in Kruskal 
coordinates until the r = 0 singularity. We include ingo¬ 
ing and outgoing velocities, and time-like, light-like and 
space-like paths. In the process we have generated the 
paths shown in |26j and more. We have determined equa¬ 
tions for the space-like caustics (degenerate paths) which 
have not been described before in this manner. We have 
used a WKB path integral approach and determined the 
Schwarzschild wave function for the collapse. Using a 
similar approach one can in the future use the paramet¬ 
ric path equations we have determined in Kruskal coor¬ 
dinates to explore the space-time inside the horizon. It 
will require expansion around the caustic. The possibil¬ 
ity of paths moving back in time might have important 
ramifications. 

It has long been argued whether quantum gravity 
should be causal or not |?7il3(j) and as of today there 
is no accepted theory of quantum gravity. We evaluate 
the propagator both for a free relativistic particle, where 
we have an exact solution, and for a collapsing ball of 
dust. Since the propagator does not vanish outside the 
light cone, we include all paths in our integration. In the 
free particle case, the WKB and exact solutions match. 
Space-like classical paths give the dominant contribution 
in the construction of the WKB approximation far out¬ 
side the light cone, where the approximation is valid. 
This propagation will be acausal (backward in time) in 
some Lorentz frames. Whether or not this proves to be 
admissible in the ultimate theory of quantum gravity is 
beyond the purpose of our work. Our work’s exploration 
of the consequence of an acausal propogator in the col¬ 
lapse of a zero pressure star might nevertheless contribute 
to the dialogue in the search of a consistent theory of 
quantum gravity. 

In addition, we propose that an understanding of these 
low probability paths that are traditionally ignored could 
shed some light in theories of black hole formation. Fur¬ 
thermore, if the black holes observed by LIGO are formed 
from scalar particles [31], the particles would be very 
light; then quantum effects might play a role in their 
formation, constraining analogue gravity theorems [321] , 
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